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Abstract 

The AMS-02 collaboration has just released the cosmic antiproton to proton ratio p/p with a 
high precision up to ~ 450 GeV. In this work, we calculate the secondary antiprotons generated 
by cosmic ray interactions with the interstellar medium taking into account the uncertainties from 
the cosmic ray propagation. The p/p ratio predicted by these processes shows some tension with 
the AMS-02 data in some regions of propagation parameters, but the excess is not signihcant. We 
then try to derive upper bounds on the dark matter annihilation cross section from the p/p data or 
signal regions favored by the data. It is shown that the constraint derived by the AMS-02 data is 
similar to that from Fermi-LAT observations of dwarf galaxies. The signal region for dark matter 
is usually required my^ 0(10) TeV and {av) ~ 0(10 cm^ s 

PACS numbers: 95.35.-|-d, 96.50.S- 
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I. INTRODUCTION 


Antimatter search in cosmic rays (CRs) are especially important for astrophysics and dark 
matter signal study. In recent years, there were great progresses in the measurements of CR 
antimatter particles. In particular, the Alpha Magnetic Spectrometer (AMS-02) launched in 
2011 has provided unprecedentedly precise measurement of CRs. Based on the results from 
AMS-02 and previous experiments, such as PAMELA and Fermi-LAT, the properties of CR 
antimatter particles have been extensively studied in a quantitative way in the literature. 

Recently, the AMS-02 collaboration has reported their measurement of the h/p ratio up to 
~ 450 GeV with a high precision [l|. The previous PAMELA p/p ratio data { 2 , 3| seemed to 
be compatible with the conventional CR prediction in the energy range of (Pfl) — (P(10) GeV, 


and can be used to set limits on the DM annihilation cross section or lifetime 


ll| . However, 


it is very interesting to note that the AMS-02 p/p ratio becomes almost flat from 10 GeV 
to several hundred GeV. This behavior is not expected from the GR secondary production 
since the measured B/G ratio decreases with energy as a power-law spectrum without any 
structure. Therefore it is necessary to perform a careful study on the secondary production 
of antiprotons taking all the uncertainties into account. 

The major challenge to interpret the AMS-02 data is how to give an accurate prediction 
of the secondary antiprotons. Although the AMS-02 collaboration has provided precise 
measurements of the proton and helium spectra, there are still large uncertainties from 
GR propagation processes. The complicated GR propagation in the Galaxy may involve 
diffusion, energy loss, convection, and reacceleration effects. The propagation parameters 
describing these effects can be determined by htting the secondary-to-primary ratios, such 


as B/G and (Sc-fTi-l-V^ 
lOBe/^Be and ^SAl/^^Al 


/Fe, and unstable-to-stable ratios of secondary particles, such as 


12 


15| . However, current measurements of these ratios are still not 


sufficient to distinguish different propagation models. Another important uncertainty comes 
from the solar modulation, which signihcantly affects low energy GR spectra, but it is not 
easy to quantify the effects. 

In this work, we give hrst a careful analysis of the secondary antiproton production. The 
GALPROP package is used to solve the GR propagation equation [ 1 ^, Q . We adopt two 
kinds of propagation models, namely the diffusion-reacceleration (DR) model and diffusion- 
convection (DG) model. Propagation parameters are determined by htting the available 
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B/C data with a Markov Chain Monte Carlo (MCMC, 17]) algorithm 18142 ll|. which is 
powerful for surveying high dimensional parameter space. Secondary antiprotons from CR 
interactions with the ISM are calculated and the uncertainties from propagation effects 
are obtained. Taking all the uncertainties into account, the secondary antiprotons seem 
marginally consistent with AMS-02 results, but show clearly different slope with energy. 

Then we study the implications of the AMS-02 p/p data on DM annihilation signals. 
Since the AMS-02 data do not unambiguously indicate an excess, we derive upper limits on 
the DM annihilation rate and signal regions where the background plus DM signals give a 
good £t to the AMS-02 data. In order to derive upper limits or signals regions, we carry out 
a global htting to the AMS-02 data with the MCMC method. Based on the AMS-02 p/p 
data, upper limits on the D M anni hilation rate and favored DM annihilation parameters 
have also been given in Refs. 


22 


25] 


This paper is organized as follows. In Sec. II, we describe the CR propagation processes. 
In Sec. Ill, we first determine the propagation parameters from the AMS-02 B/C ratio data 
and then calculate the secondary antiproton flux. In Sec. IV, we study implication on DM 
annihilation by the AMS-02 p/p ratio data. Finally, we give the discussions and conclusions 
in Sec. V. 

II. PROPAGATION OF GALACTIC COSMIC RAYS 


After accelerated in sources. Galactic CRs are injected and diffuse in the interstellar space 
and suffer from several propagation effects before they arrive at the Earth. A conventional 
assumption is that CRs propagate in a cylindrical halo with a half height Zh, beyond which 


CRs escape freely. The propagation equation can be expressed as 


26] 


^ = Q(x,p) + V ■ - Ve^) + ^ 


' 2n 

P ^PPdpp2 


dp 


Pp! 




t. 


p) 


Tr 


( 1 ) 


where p) = pj{'x,p,t) is the CR density per momentum interval, Q(x, p) is the source term, 
p = dp/dt is the momentum loss rate, and the time scales r/ and characterize fragmen¬ 
tation processes and radioactive decays, respectively. 

The convection velocity is usually assumed to linearly depend on the distance away 
from the Galactic disk and the convection effect can be described by the quantity dVc/dz. 
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TABLE I. Propagation models considered in this paper. 


Model Propagation parameters 


Description 


DR Dq, 5, VA, Zh Diffusion + reacceleration, r] = 

DR-2 Dq, S, va, Zh Diffusion -|- reacceleration, rj = - 

DC Dq, 6, dVc/dz, Rq, Zh Diffusion -|- convection, rj = 1, 6 


The spatial diffusion coefficient D^x can be parameterized as [2 


1, Ro = 4 GV 
0.4, Ro = 4 GV 
= 0 for R < Ro 


Dxcc = DqI3^ {R/RqY , 


( 2 ) 


where R = pc/Ze is the rigidity, (3 is the CR particle velocity in units of the light speed c, and 
Dq is a normalization parameter. Although the slop of diffusion coefficient 6 is predicted to 
be (5 = 1/3 for a Kolmogorov spectrum of interstellar turbulence, or 5 = 1/2 for a Kraichnan 
cascade, it is usually taken as a free parameter when explaining data. The factor of 
denotes the effect that the diffusion coefficient could be altered at low velocities due to the 
turbulence dissipation. 

The CR reacceleration process due to collisions with interstellar random weak hydrody¬ 
namic waves can be described by the diffusion in momentum space with a coefficient Dpp, 
which is related with Dxx by 281 


DppDxx 




h2)(4-h)n;’ 


( 3 ) 


where va is the Alfven velocity and u: is the ratio of the magnetohydrodynamic wave energy 
density to the magnetic held energy density. Since Dpp is proportional to v\/u, we can 
set a; = 1 and just use va to characterize the reacceleration effect. Therefore, the major 
propagation parameters involve Dq, 5, Rq, p, va, dVc/dz, and Zh- 

When CRs propagate in the solar system, their spectra with R < 20 GV are signih- 


cantly affected by the solar wind. This is the solar modulation effect, which depends on 
the solar activity and varies with the solar cycle. It can be described by the force held 


approximation 


29l | with only one parameter, i.e., the solar modulation potential 0. 


Due to fragmentation and radioactive decays, secondary CR particles are produced in 
propagation processes. As a result, secondary-to-primary ratios and unstable-to-stable sec¬ 
ondary ratios, such as B/C and ^°Be/®Be, are sensitive to propagation paramors, but almost 
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independent of primary injection spectra. Thus the measurements of these ratios are useful 


to determine the propagation parameters (see e.g. Refs. jSOl. l31|). 


Note that there are degeneracies between the propagation models with reacceleration 
process and those with convection effects. Current measurements are not sufficient to dis¬ 
tinguish these models. Therefore, we separately consider the DR and DC models when 
calculating the antiproton flux. The descriptions of these models are given in Table HI Ro is 
taken to be 4 GV in the DR and DR-2 models. In the DC model, however, Ro is set to be 
a free parameter and 5 = 0 is imposed for R < Rq. 

III. ASTROPHYSICAL PREDICTION FOR THE p/p RATIO 


In this section, we calculate the secondary antiproton flux generated by CRs when they 
propagate in the Galaxy considering the uncertainties from propagation processes. In the 
following, we first derive the ranges of propagation parameters. 

It is generally assumed that GR particles are accelerated in supernova remnants (SNRs). 
Thus we assume the spatial distribution of the GR sources follow the SNR distribution as 

b{r — r©) 


f{r,z) = 


re 


exp 


rQ 


exp — 


( 4 ) 


where the distance from the Sun to the Galactic Genter Tq = 8.5 kpc, and the characteristic 


height of the Galactic disk Zg ~ 0.2 kpc. We adopt a = 1.25 and b = 3.56 


31 | in the 


calculation. In the DR-2 and DG models, the injection spectra of primary nuclei are taken 
to be single power-law functions. In the DR model, the injection spectra are assumed to be 
broken power-law functions with respect to the rigidity as 

-Ul 



/ R 

qi oc 1 

TT- 

\Rhv 


( ^ 

qi oc 1 

-w- 

V-^br 


- 1/2 


( 5 ) 


, R > Rbr- 


Then the injection source term is Qi{'x,p) = f (r,z)Qi(p). 

We use the numerical tool GALPROP 13|, ll6j to solve the propagation equation. The 
propagation parameters are determined by fitting to the nuclear data B/G and ^°Be/®Be. In 
order to improve the htting efficiency, the MGMG method is employed to derive the posterior 
probability distributions of the propagation parameters. By fitting the B/G ratio data from 


AMS-02 


32 | and AGE 33|, and the ^°Be/®Be ratio data from AGE 34], Balloon 35-37]. 
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TABLE II. Mean values and la uncertainties of the propagation parameters derived through fitting 
the data of B/C and ^*^Be/®Be ratios in three propagation models. 



DR 

DR-2 

DC 

Dq (10^® cm^ s 

6.58 ±1.27 

3.59 ±0.88 

1.95 ±0.50 

6 

0.333 ±0.011 

0.423 ±0.017 

0.510 ±0.034 

Ro (GV) 

4 

4 

4.71 ±0.80 

VA (km s“^) 

37.8 ±2.7 

22.6 ±3.1 

/ 

dVc/dz (km s“^ kpc“^) 

/ 

/ 

4.2 ±3.2 

Zh (kpc) 

4.7 ± 1.0 

3.5 ±0.8 

2.5 ±0.7 

0B/C (MV) 

326 ± 36 

334 ± 37 

182 ± 25 


IMP7&8 


38|, ISEE3-HKH |39|, ISOMAX |40|, Ulysses-HET |4l|, and Voyager [421, we 


obtain the mean valnes and la errors of the propagation parameters for the three models, 
as shown in Table m We also show the fitting resnlts of B/C ratio within 95% confidence 
ranges compared with the experimental data in Fig. [H 



FIG. 1. Fitting results of the B/C ratio within 95% confidence ranges in three propagation models, 


compared with the data from AMS-02 


32 1 and ACE 


33|. 


Then the spectrum of secondary antiprotons and p/p ratio are calculated taking into 
account the propagation parameter uncertainties. Here the proton injection spectrum is 
determined by htting to the AMS-02 321 and CREAM data. It should be noted that 
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the proton spectrum hardening at ~ 200 GeV enhances the p/p ratio at the high energy end. 
But even with such proton spectrum hardening the predicted p/p ratio shows a trend that 
decrease with energy for E > 100 GeV, which seems not consistent with the AMS-02 data. 
The p/p ratio predictions for propagation parameters varying within 95% conhdence ranges 
in the three propagation models are presented in Fig. [21 The solar modulation potential 0 
are taken as 0.7 GV for the DR and DR-2 models, and 0.5 GV for the DG model. 


We notice that the DR model underproduces antiprotons at low energies compared with 
the AMS-02 data. Actually such an underproduction has been realized for a long time when 


studying the PAMELA 


, y] and BESS 44] results 


9 


31 


45[ |. In order to solve this 


problem, the diffusion coefficient at low energies should be modified 45j. Other effects, such 
as the charge-dependent solar modulation and convection, can also be used to avoid this 


discrepancy j^, [ll|. From Fig. [21 we can see that the DR-2 model with rj = —0.4 and the 
DG model with a break in the diffusion coefficient can well ht the data up to ~ 100 GeV. 
For Ek > 100 GeV, however, the predicted p/p ratios in all these models tend to decrease, 
and seem to be hardly reconciled with the AMS-02 data. 



FIG. 2. Astrophysical predictions of the p/p ratio for propagation parameters varying within 95% 
confidence ranges in the DR (red band), DR-2 (blue band), and DC (green band) propagation 
models. Data points denote the AMS-02 p/p ratio measurement [l|. 
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IV. IMPLICATIONS FOR DARK MATTER ANNIHILATIONS 


In this section we stndy the implications of the AMS-02 p/p data on DM annihilation. 
As the AMS-02 data at high energies and the propagation model uncertainties are still 
large, there is no obvious antiproton excess in data. However, the high energy discrepancies 
between the predicted astrophysical background and the data leave space for contributions 
from extra antiproton sources. It is usually expected that DM particles in the Galactic 
halo may annihilate into quarks and gauge bosons, and produce antiprotons after decay and 
hadronization processes. 

For self-conjugate DM particles, the source term of antiprotons from DM annihilations 
is given by 




2 


( 6 ) 


where is the DM particle mass, and (cm) is the thermally averaged annihilation cross 
section. dN/dE^ is the differential number of antiprotons produced per annihilation. We 
use the results of PPPC 4 DM ID 461 including the electroweak corrections to calculate 
dN/dEk- The DM density distribution p(r) is adopted to be the Navarro-Frenk-White 
(NFW) prohle Idsj 


p{r) = 


Ps 


(7) 


(r/rj(l -F r/rsY' 

where = 20 kpc, and ps is derived from the condition that the local DM density at 8.5 kpc 
from the Galactic Genter is 0.3 GeV cm“^. 

When studying the DM implications we take the best Pi p/p background from the bands 
in Fig. |2]and then add the contributions from DM annihilations to the p/p ratio. Since the 
p/p ratio prediction in the DR model seems inconsistent to the data, we only consider the 
DR-2 and DG models below. In order to give the best £t to the AMS-02 data we introduce 
a free factor Cp to normalize the flux of secondary antiprotons. This factor represents pos¬ 
sible uncertainties in astrophysical antiproton production processes, such as uncertainties 


from hadronic interactions (see e.g. Refs. l50|), ISM density distributions, and nuclear 


enhancement factors from heavy elements. The solar modulation potential (pp is another 
free parameter to give best £t to data. 

Then we derive the upper limits on the DM annihilation cross section at 95% G.L. from 
the AMS-02 data taking the uncertainties from propagation into account, as shown in Fig. 
[3l Here we consider DM annihilations to bb and W~^W~ hnal states, which are typical 
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cases for antiproton prodnctions. The snm of the contribntions from backgronnd and DM 
annihilation to the antiproton flnx is given by 


(crn), Cp, 0p) = <hbkg(cp, 4>p) + {crv ), 0p). 


( 8 ) 


For each valne of we derive an npper limit on {av) by reqniring 


^^((cTn)) -xLn = 3.84, 


(9) 


where Xmm corresponds to the best £t. 






FIG. 3. The contours denote the 68% and 95% confidence regions favored by the AMS-02 p/p data 
in the my.-{crv) plane for the bb (upper panels) and W~^W~ (lower panels) annihilation channels. 
The red band represent the upper limits at 95% C.L. on the DM annihilation cross section set 
by the AMS-02 p/p data taking the uncertainties from propagation. For a comparison, the upper 
limits from the Fermi-LAT observations of dwarf galaxies are also shown as blue solid lines. The 
left and right panels correspond to the results in the DR-2 and DC models, respectively. 
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TABLE III. Fitting results for DM annihilations into bb and W~^W in the DR-2 and DC propa¬ 


gation models. 


bb channel 

log(m^/GeV) 

log(((T?;) /cia? s 

Cp 

h (GV) 

X^/dof 

DR-2 

best 

mean 

4.65 

4.29 ± 0.45 

-21.97 

-22.57 ±0.72 

0.930 

0.928 ±0.013 

0.723 

0.724 ± 0.030 

11.8/26 

DC 

best 

mean 

4.99 

4.48 ± 0.37 

-21.02 

-21.89 ±0.61 

0.945 

0.944 ±0.012 

0.696 

0.696 ± 0.030 

17.6/26 

w+w- 

channel 

log(m^/GeV) 

\og{{(Tv) /CW? S“^) 

Cp 

h (GV) 

X^/dof 

DR-2 

best 

mean 

3.85 

4.10 ±0.46 

-23.05 

-22.62 ±0.81 

0.932 

0.934 ±0.011 

0.719 

0.717 ±0.029 

10.90/26 

DC 

best 

mean 

3.94 

4.15 ±0.38 

-22.50 

-22.15 ±0.69 

0.947 

0.947 ±0.011 

0.691 

0.691 ±0.029 

15.7/26 


From Fig. [3] we can see that the upper limits in bb and W~^W~ channels are very similar 
due to their analogous initial antiproton spectra. The DR-2 propagation model gives slightly 
stricter constraints, compared to the DC model. 

Since there are discrepancies between the predicted astrophysical backgrounds and the 
data at high energies, DM contributions to the antiproton flux can signiflcantly improve the 
fitting to the AMS-02 p/p data. We carry out a global fitting by varying {cp, (j)p, m^, (an)}, 
in which is set to below 100 TeV, in order to avoid contradicting with the unitarity 
bound. Then we obtain the best fit values, means values, and la uncertainties of these 
parameters, as listed in Table IIIII The p/p ratio for the best-fit DM parameters are shown 
in Fig. m We can see that ~ (9(10) TeV and (an) ~ (9(10“^^) cm^ s“^ could explain 
the data very well. 

Fig. |3]also shows the favored regions at 68% and 95% C.L. in the m^-(an) plane. The 
favored values of (an) are highly correlated with the values of m^, since the DM source 
term is proportional to (an) /m^, as indicated in Eq. ([6]). Since the DM annihilations into 
quarks or gauge bosons also produce high energy photons, y-ray observations, such as those 
from Fermi-LAT, would set strong constraints on the DM annihilation cross section. Also 
shown in Fig. |3] are the upper limits from the observations of spheroidal dwarf galaxies with 


six-year Fermi-LAT gamma-ray data 


5l| . It can be found that the constraint from AMS-02 
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FIG. 4. p/p ratios including the contributions from DM annihilations into bb (red lines) and W~^W 
(blue lines) for our best-fit parameters. Data points denote the AMS-02 p/p ratio measurement 
The left (right) panel corresponds to the results in the DR-2 (DC) model. 


B. 


is similar to that derived from Fermi-LAT dwarf galaxy observation. 


V. CONCLUSIONS 

In this work we investigate the new AMS-02 resnlts of p/p ratio. In order to nnderstand 
the data, the crucial step is to perform a reliable study on the secondary antiproton flux 
generated by CR interactions with the ISM. We consider three propagation models, namely 
the DR, DR-2 and DC models dehned in Tab. [B The prorogation parameters are constrained 
by the B/C data using a MCMC algorithm. We hnd that all the models can provide very 
good £t to the B/C data. 

Within la ranges of the fitted propagation parameters, we give the predictions of the 
secondary antiprotons in the three propagation models and compared them with the AMS- 
02 p/p data. It is found that the DR model underproduces antiprotons at low energies, 
which has already been noticed for a long time. The DR-2 and DC models can give a quite 
good £t to the AMS-02 data below ~ 100 GeV. However, it should be emphasized that the 
predicted secondary p/p ratio above ~ 100 GeV has a different slop from the AMS-02 data. 
The data keep flat up to ~ 450 GeV, but the predictions of all the three models decrease 
with the energy. The key factor is that the diffusion coefficient is in a form of single power 
law, namely D^x ~ R~^, which is determined by the B/G data. If the diffusion coefficient 
Dxx permits a ffatter behavior at high energies within the allowance of the B/G data, a 
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flatter behavior of p/p could also be expected [52 1. 

If we consider an extra source of antiprotons from the Galactic DM annihilation, the 
htting to the AMS-02 data can be signihcantly improved. The data can be well explained by 
a DM particle with a mass of ~ lOTeV and an annihilation cross section of ~ 10“^^ cm^ s“^ 
into the bb or W~^W~ hnal states. 

If we assume that there is no obvious excess in the p/p ratio, the AMS-02 data can be 
used to set upper limits on the DM annihilation cross section into bb or W~^W~. We hnd 
that the constraints given by the DR-2 model are slightly stricter than those given by the DC 
model. For small m^, the constraints from AMS-02 are even stronger than the constraints 
from Fermi-LAT. 
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